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1 Introduction 



The B-meson distribution amplitude was introduced in P] as the direct analogue of light- 
cone distribution amplitudes of light mesons [3 El El in an attempt to describe generic 
exclusive B-decays by the contribution of the hard gluon exchange. Since then, consid- 
erable effort has been invested to understand the QCD dynamics of heavy meson decays 
in the heavy quark limit. The radiative decay B ^eu provides one with the simplest 
example of such processes jS] . This form factor can be calculated in terms of the B-meson 
distribution amplitude to one-loop accuracy [Oj and arguments have been given that the 
corresponding factorization formula is valid to all orders in the strong coupling 0. Similar 
QCD factorization formulas have also been proposed for the related processes S — 77 
and B —>■ '~fi~^i~ |H1 E]. For weak decays involving energetic light hadrons in the final 
states the QCD factorization is more complex since one must isolate the end-point soft 
contributions in terms of additive contributions. This is a hot topic, see e.g. ^U]) and the 
results have been encouraging although, as it has been repeatedly pointed out [TTl [T2] , 
the l/rrib corrections to heavy-to-light exclusive decays are most likely large and require 
quantitative treatment. 

The B-meson distribution amplitude plays the central role in all known factorization 
formulas, but, surprisingly, received relatively little attention in the past. In the present 
work we use QCD sum rules to present a realistic model for the B-meson distribution 
amplitude, consistent with all QCD constraints. We also take this opportunity to clarify 
its theoretical status on which there has been certain confusion. The approach that we take 
in this paper is inspired by the classical work by Chernyak and Zhitnitsky on the QCD sum 
rule analysis of the distribution amplitudes of light mesons [13]. In particular we argue 
that the relevant matrix element of a bilocal quark-antiquark operator can be calculated 
by the QCD-corrected expansion at imaginary light-cone distances. We find that the 
nonperturbative corrections remain under control and present quantitative estimates for 
the distribution amplitude and its first inverse moment which enters decay form factors 
at tree level. Our results can be considered as an extension of an earlier QCD sum rule 
calculation by Grozin and Neubert [21 (see also with the main difference being that 
we calculate NLO radiative corrections to the sum rule. This is important since the true 
analytic structure of the distribution amplitude is only seen at this level. 

The presentation is organized as follows. Sec. 2 is introductory and contains the nec- 
essary definitions. A simple model for the distribution is obtained in Sec. 3 using QCD 
perturbation theory and duality. The complete sum rule is constructed in Sec. 4 where we 
discuss the structure of nonperturbative corrections. This Section also contains our main 
results, the summary and conclusions. 

2 Definitions 

Following Ref. ^B] we define the B-meson distribution amplitude as the renormalized 
matrix element of the bilocal operator built of an effective heavy quark field hy{0) and a 
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light antiquark q{tn) at a light-like separation: 



(0| q{tn)fl[tn,0]Th,,{0) \B{v)) 



R 



with 



[tn, 0] = Pexp 



(2) 



Here is the heavy quark velocity, is the light-like vector, = 0, such that n - v = 1, 
P+ = is the projector on upper components of the heavy quark spinor, T stands 

for an arbitrary Dirac structure, \B{v)) is the 5-meson state in the heavy quark effective 
theory (HQET) and is the decay constant in HQET, which is related to the physical 
B-meson decay constant to one-loop accuracy as 



fBy/rriB = F{^) 
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The notation [. . .]r in stands for the renormalization in a MS-like scheme and /z here 
and below refers to the MS normalization scale. 

The invariant function where t is a real number defines what is usually called 

the leading twist B-meson distribution amplitude in position space, in contrast to the 
amplitude /i) which involves a different light-cone projector — ^ instead of ^ — in 
between the quarks; here = 0,n ■ n = 2. This name is not exact since the translation 
symmetry of the theory is broken by presence of the effective heavy quark field and hence 
neither geometrical nor coUinear twist are defined. In the present paper we only consider 
the distribution ^j^{t,^) and its Fourier transform 
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where in the first equation the integration contour goes below the singularities of {t, fi) 
that are located in the upper-half plane. 

The scale dependence of the distribution amplitude is driven by the renormalization 
of the corresponding nonlocal operator 0+(t) = q(tn) 'i/[tn,0]T h^{0). The corresponding 
Z-factor was computed in p!F to one-loop order. In momentum space, the result reads 
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where 
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Figure 1: One-loop renormalization of the nonlocal light-cone operator built of one light and 
one effective heavy quark field (double line). The dashed line indicates the gluon Wilson line 
insertion in between the quark fields. 



with d = 4 — e and the standard notation 2/e = 2/e — 7£; + ln47r. Here [...]+ is the usual 
"plus" -distribution. 

In order to understand the meaning of this result it is instructive to consider operator 
renormalization in position space. The corresponding to (jU)) one-loop expression is 

0;-(t,/i) = Ol-^it) + ^ { + JH^t^^)^ Ol-^it) 

€ Jq 1 — u e 

where the first two terms in curly brackets correspond to vertex-type corrections shown 
in Fig. la and Fig. lb, respectively, (in Feynman gauge) and the third term takes into 
account the quark field renormalization: q'^^^ = Zq^^'^q"^'^^ and = Z^^^'^h}^^^'^ with 
Zq = l- (2/e)(a^CF/47r), Zh = l + {A/e){asCF/A'K) |T7|. The exchange diagram in Fig. Ic 
is UV- finite and does not contribute [HI HH] . 

Note the following property: renormalization of the nonlocal light-cone operator 0+(t) 
((Tj) is quasilocal: it only gets mixed with itself and with operators with smaller light- 
cone separation. In fact the heavy quark vertex correction in Fig. la corresponds to a 
multiplicative (cusp) renormalization in coordinate space ^M, while the light quark vertex 
correction is identical to the similar contribution to the light-quark-antiquark nonlocal 
operators [12]. For light quarks, this property of quasilocality guarantees existence of the 
Wilson short distance operator product expansion (OPE) since it implies that the "size" 
of the operator is not altered by the renormalization. In the present case, however, the 
local OPE does not exist because of the term \n{it^) which is non-analytic at t — > 0. It 
is easy to see that this contribution arizes from the term ~ {fity/e'^ in the dimensionally 
regularized diagram in Fig. la so that the answer for this diagram depends on the order 
of limits t and 4 — = e — * 0. We conclude that renormalization of the nonlocal 
light-cone operator built of one light and one effective heavy quark field does not commute 
with the short- distance expansion. In particular 

[qitn)ji[tn,0]TKiO)]R ^ $^-[^(0)(Z) -nYKmR, (8) 

p=0 ^' 
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and the equality does not hold even as an asymptotic expansion. As a consequence, non- 
negative moments of the B-meson distribution amplitude J dk kP(j)^{k, fi) for p = 0, 1, 2, . . . 
are not related to matrix elements of local operators and in fact do not exist: It is easy 
to see that the logarithmic singularity of the amplitude in position space /u) ~ In (it) 
for t ^ implies that the Fourier integral (jH) is logarithmically divergent at /c — oo, 
that is (j)^{k,fi) ~ for ^ /i, in agreement with the analysis in ^H]- The analysis 
of moments Jdk k'^(j)^{k, fi) in [20] tacitly assumes a different definition of the B- 
meson distribution amplitude, such that the nonlocal light-cone operator is defined as the 
generating function for renormalized local operators on the r.h.s. of (jH])- This implies 
e.g. that power divergences are subtracted. This is a different object which, most likely, 
does not satisfy any simple renormalization group equation and has no obvious relation 
to exclusive B-decays. 



3 Sum Rules: Perturbation Theory 

Aim of the present study is to suggest a realistic model of the B-meson distribution am- 
plitude that would be consistent with all QCD constraints. To this end we evaluate the 
necessary B-meson matrix elements using the standard QCD sum rule approach [21]. In 
this section we set up the framework and present intermediate results that only include 
perturbation theory contributions and the assumption of duality. The complete treatment 
including nonperturbative corrections is presented in the next section. 

To derive the sum rules we consider the following correlation functions in HQET: 

I j d'x e-''^^^^\0\T{m^M0)h.{x)T2q{x)}\0) = -^Ti [T,P+T 2] U{lj) (9) 

and 

I J d^x e-'^^^^'^OlTiqitn) ]iT^[tn,0]KiO)K{x)T2qix)}\0) = 

(10) 

The correlation function Il{uj) has a pole at u = A where A = niB — mi, is the usual HQET 
parameter, and the residue at this pole is proportional to the HQET decay constant -F(/i): 

H(co') = -F^(/i) h higher resonances and continuum. (11) 

2 A — a; 

Similarly, 

Tit, u) = -F\fi) J / dk e-'^* Mk, /x) + • • • (12) 

2 A-u Jo 

On the other hand, both correlation functions can be calculated in QCD at negative values 
of u of the order of 1 GeV in perturbation theory and taking into account nonperturba- 
tive effects induced by vacuum condensates [21]. Matching the two representations one 
obtains a sum rule. There are two technical details: First, one makes an assumption that 
contributions of the continuum and of higher resonances can be taken into account by the 



4 



restriction to the so-called duality region < s < ujq in the dispersion representation for 
the correlation functions, e.g. 



ds 



UJ 



■Pu[s) 
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where pn{s) is the corresponding spectral density. The numerical value for the parameter 
ujq (called continuum threshold) is usually taken to be in the interval 0.8 — 1.0 GeV 
|22[ 1^ 1^ II 7j. Second, one makes the so-called Borel transformation 



/ -^Pnis)^ / ds e-^/'' Puis) 
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introducing the variable M (Borel parameter) instead of the energy uj in order to suppress 
higher-order nonperturbative corrections and minimize the dependence on the continuum 
model. The resulting sum rule for the correlation function Il{w) is well known 
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Here = ctsip), {qq) — —(240 MeV)'^ is the quark condensate and ml is the ratio of the 
quark-gluon and quark condensates ml = {qg{crG)q) / {qq) ~ 0.8 GeV^. The given numbers 
correspond to the renormalization scale /i = 1 GeV. With the choice A = 0.4 — 0.5 GeV 
and Wo = 0.8 — 1.0 GeV the sum rule in (fTHjl is satisfied for a wide range of values of 
the Borel parameter 0.3 GeV< M < oo and is used j221 120] to determine the B-meson 
decay constant F{fi) in the heavy quark limit. In the numerical estimates in this paper we 
will take the "window" 0.3 GeV < M < 0.6 GeV in which the matching is done [21] and 
use the value ^^(l GeV) = 0.5 (Aq^'^^^ ~ 360 MeV) which is consistent with the world 
average. 

Our task in this work is to derive the similar sum rule for the correlation function T{t, u) 
defined in fllOp . The perturbative contributions are shown in Fig. 2. The corresponding 
sum rule reads, so far without nonperturbative corrections: 
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Figure 2: Correlation function in QCD perturbation theory to first order. 

p>(s, k, li) = ^'^^J^^ \-x + ln(l -x)- 2(1 - x) ln(l - x) + 2 \n^{l - x) 
A; 

+21n-[x + ln(l -x)] 
A* 



(17) 



Here Li2(x) is Euler dilogarithm function and we used a shorthand notation x = 2s/ k, 
Neglecting corrections for a moment, one gets a simple expression 
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lF'ifi)e-~^/'' 0+(A;,/x) = ^9{2wo - k) k I ds e"^/^ . [L^) 

In the local duality limit M —>■ oo using the sum rule expression for F{n) ()15|) with the 
same accuracy, l/2F^(yu) ~ {Nc/67!''^)u!q, one obtains 
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(jyJkr'' = — ^^(2^0 - k) k(2uJo - k) (19) 

which reminds the asymptotic light-cone distribution amplitude of light mesons if rewritten 
in terms of the scaling variable ^ = k/{2ujo)- For finite values of the Borel parameter M 
the B-meson distribution amplitude gets skewed towards smaller values of the momentum 
but qualitatively remains the same, see Fig. 3. Note that it has finite support k < 2ujo 
and can be interpreted as the probability amplitude to find the light quark (on-shell) in 
the B-meson with momentum k. 

Beyond the Born approximation this simple parton-model interpretation is lost since 
the distribution amplitude develops a high-momentum "tail" with k > 2uJ() and in this 
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Figure 3: B-meson distribution amplitude (/>_|_(A;,// = 1 GeV) calculated from the sum rule 
in QCD perturbation theory to leading order (dashed curves) and next-to-leading order 
(solid curves) for the continuum threshold ujo=1 GeV and two values of the Borel parameter 
M = 0.3 GeV (left panel) and M = 0.6 GeV (right panel). The value of the decay constant -F(/u) 
appearing on the l.h.s. of the sum rule is substituted by the corresponding sum rule (|15|) with 
the appropriate accuracy (LO or NLO) and neglecting the condensate contributions. 

region cannot be thought of as a probability amplitude for the two-particle state on mass 
shell. The 0{as) radiative correction turns out to be very large (~100% of the Born term) 
but cancels to a large extent against the similar large radiative correction to F{n) pK|l22j. 

The numerical results for two values of the Borel parameter M = 0.3 GeV and M = 
0.6 GeV are shown in Fig. 3. We choose i^o = 1 GeV for this plot and substitute the 
coupling -F^(/u) appearing on the l.h.s. of the sum rule ()16|) by the sum rule (fT3j) to the 
same accuracy, i.e. neglecting nonperturbative corrections. In this way the dependence 
on A cancels out and the sensitivity to other parameters {ujq and M) is strongly reduced. 
Indeed, it is seen in Fig. 3 that dependence on the Borel parameter is rather mild. Note 
that for large k the distribution amplitude becomes negative. The asymptotic behavior is 

(f)Jk) r^k ioT k^O, ^+(k) — ln(A://i) for > /i , (20) 

k 

in agreement with [16 . Also the scale dependence of the distribution amplitude extracted 
from the sum rule (fT^ agrees with ^Sl- All results are shown for /i = 1 GeV. 

Of particular interest for the QCD description of B-decays is the value of the first 
negative moment 

lik 

(21) 

We obtain from the sum rules: 

A-i = 1.49 - 1.83 GeV"^ for cuq = 1.0 GeV, 

= 1.79 - 2.08 GeV"^ for tuo = 0.8 GeV , (22) 

where the lower value corresponds to M = 0.6 GeV and the higher one to M = 0.3 GeV 
for each choice of ujq. Notice that A^^ decreases as M increases and in the local duality 
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limit we obtain 

where it is taken into account that in the hmit M ^ oo the sum rules effectively become 
normalized at the scale /i = 2ujq because of subtraction of the continuum from the running 
coupling, cf. |.26J. To avoid misunderstanding we remind that all results of this section 
correspond to the sum rules in QCD perturbation theory and the given numbers will be 
superseded by those in the next section where we consider the nonperturbative corrections. 
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(23) 



4 Sum Rules: Nonperturbative Corrections 

The primary source of nonperturbative corrections to the sum rules in HQET is provided 
by the quark condensate. The corresponding diagrams (leading and next-to-leading order) 
are shown in Fig. 4. The leading-order contribution in Fig. 4a is simply 

T(^i)^t,uj) = ^. (24) 
2co' 

It does not depend on the quark- ant iquark separation and gives rize to the 5-function type 
contribution to the r.h.s. of the sum rule in (fT^ 

. . . - i {qq)6{k) . (25) 

Condensates of higher dimension produce even more singular terms, the expansion goes in 
derivatives of the (5-function at /c = 0. This is a well-known problem which is familiar from 
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the QCD sum rule studies of light-cone distribution amplitudes of light mesons [131 ^^^1 
nucleon parton distributions ^71 |2H] '■ The short- distance OPE which is the the basis of 
the SVZ approach is inadequate for a calculation of distribution functions point-by-point 
in the momentum fraction space. As a consequence, QCD sum rules cannot be used for a 
direct calculation of distribution amplitudes (unless they are supplemented by additional 
assumptions) but rather provide constraints which have to be implemented within QCD- 
motivated parameterizations (models) of the distributions, consistent with perturbative 
QCD. For a model to be self-consistent, there are three conditions: 

• The end-point behavior of the distributions has to be consistent with QCD. 

• The model has to be closed under the QCD evolution, i.e. calculation of the scale 
dependence has to be possible and not involve nonperturbative parameters other 
than those specified by the model at the reference scale. 

• The model has to involve a minimum possible number of nonperturbative parame- 
ters. 

The Chernyak-Zhitnitsky models of light-cone distribution amplitudes of light mesons 
give the classical example of such an approach. In this case one expands the distribution 
amplitude in a series over orthogonal polynomials, e.g. for the pion 

oo 

0^(x,/i)=6x(l-x) Mf^)C'J\2x-l), (26) 

p=0,2... 

SO that coefficients in this expansion correspond to (Gegenbauer) moments of (pT^ix), and 
defines a model by truncating this expansion at a certain p = Pmax- The first pmax coef- 
ficients are then estimated using QCD sum rules. (In practice one takes pmax = 2 since 
estimates of higher-order coefficients turn out to be unreliable.) The model satisfies all the 
above criteria since the correct end-point behavior is built in by construction and higher- 
order coefficients can only get mixed with lower-order coefficients but not vice versa; it 
follows that the set of coefficients {ipo, ip2, ■ ■ ■ , V'fcmax} closed under renormalization''' and 
the distribution amplitude (j).j,{x,fi) can be calculated at arbitrary scale from its model at 
IJ, = Ho- It indeed involves a minimum number of parameters, each of which has a clear 
meaning in QCD as the matrix element of a certain local operator and can eventually be 
calculated e.g. on the lattice. 

In contrast to (j26|) . the B-meson distribution amplitude cannot be written as a sum 
of independent terms that have autonomous QCD evolution but rather is given by the 
integral in the complex moments plane |16^. This feature reminds evolution of parton 
distributions in the deep-inelastic inclusive lepton-hadron scattering, but in difference to 
the latter case one cannot obtain complex moments of the B-meson distribution amplitude 
by analytic continuation from the set of real integers (as we mentioned in Sec. 2, every 
non-negative moment of (p+{k, /i) diverges). As the result, one necessarily has a continuous 
rather than discrete set of nonperturbative parameters. 

^The coefficients in the Gegenbauer expansion are renormaUzed multipficatively to leading order be- 
cause of conformal symmetry; tliis property is, however, not essential for our argument. 
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One option |2Hl 1211 Ull is to parametrize the B-meson distribution amplitude by the 
matrix element of the bilocal operator in at imaginary light-cone separation 

t = -ir, ^+{t, fi) = <^+{-iT, fi) . (27) 

Obviously 

ip4r,f^)= / dke'^''Mk,fi) (28) 
Jo 

and the parameter As is given by the simple integral 

POO 

A^i(/i)= / dr<^+(r,/i). (29) 



The purpose of going over to imaginary light-cone times (distances) is similar to that of the 
usual Wick rotation: In this way the oscillating exponents corresponding to the light-cone 
time dependence of intermediate states propagating along the light-cone are converted to 
falling exponents suppressed by the energy of the state, where the light-cone quantization 
is implied. Simultaneously, the normalization scale fi acquires the physical meaning of 
the cutoff in energy of the intermediate states. Note that the renormalization of v9+(r, /x) 
only involves the distribution at smaller light-cone separations, cf. (|7j). This implies that 
knowledge of (y9+(r, /i) at small distances up to r < Tmax is sufficient to calculate its scale 
dependence in the same distance range, in agreement with the self-consistency criterium 
formulated above. 

On the other hand, it is easy to understand that the function {p+{T,fi) can be calcu- 
lated at small r using OPE; expansion in vacuum condensates of increasing dimension 
corresponds to the Laurent expansion of </9_|_(r, /i) in powers of r, which is modified by 
calculable perturbative corrections. The condensate expansion seems to be under control 
up to distances of order r ~ 1 GeV^^ (that is, of order 0.2 fm). At larger distances the 
OPE diverges and one has to either truncate f+{T, /i) at a certain Tmax or rely on a certain 
model for the large r behaviour. Provided that the nonperturbative corrections decrease 
sufficiently fast for large r one can hope that the model assumptions do not lead to a large 
uncertainty in the overall result. 

To illustrate this construction, we have calculated the quark condensate contribution 
including the a^-correction, see Fig. 4, the contribution of the gluon condensate. Fig. 5, 
and the contribution of the mixed condensate (qcrgGq) ~ ml{qq) which is obtained as the 
expansion of the diagram in Fig. 4a in the background gluon field. The resulting sum rule 
in which we have also included the perturbative contribution of Fig. 2 reads 

V(/.)e-^/^^^+(r,/i) = (30) 
da e-/'"pp„.(s,r,,,) - -(qq)\l + - — - W(Tpe™) - ln(T^e™) 

where the perturbative spectral density can be read off Eq. The contribution of the 

gluon condensate turns out to be very small and will be neglected in what follows. We 
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tn 

Figure 5: Gluon condensate contribution to the correlation function (|lfl|) . Only this diagram 
contributes in the Fock-Schwinger gauge. 

further note that the sum of the diagram in Fig. 4d and one half of the heavy quark self- 
energy correction in Fig. 4f define the universal renormalization factor of the Wilson line 
built of the light-like segment of length —ir and the time-like segment of length —i/M. It 
can be shown that the corresponding contributions ~ (a^Ci?)" exponentiate to all orders 
PU] and produce a Sudakov-like exponential suppression factor 



S{t, M, fi) = exp 



27r 



In^ (rue^^) + — - 1 - In + L2(-2rM) 

^ ^ ^ 24 2M ' 



(31) 



which is the same for the quark and the quark-gluon condensate. Note that L2(— 2rM) ~ 
— ^ln^(2rM) for r ^ 1/M. We end up with an improved sum rule 

^F2(/i)e-^/^V^.(r,/i) = (32) 
= j^s e-/^Vpert(s,r,/x) - \{qq) S{T,M,fi) |l + ^[2 - ln(r/.e^-) " In ^ 

-ln(l + 2rM)]--^(l + 2rM) 

in which the double-logarithmic corrections to the quark and quark-gluon chiral conden- 
sates are resummed. We do not attempt the similar resummation in the perturbative 
contribution since its effect is negligible compared to the 1/t^ falloff inherited from the 
Born term. 

The perturbative and the nonperturbative contributions to the sum rule result ()32|) for 
V2+(t, fi) are shown separately as a function of distance r in Fig. 6 (cuq = 1 GeV) and Fig. 7 
{ujq = 0.8 GeV) for two different values of the Borel parameter M = 0.3 GeV and M = 
0.6 GeV. Note that at small distances the nonperturbative corrections are significantly 
smaller than the perturbative contribution. The nonperturbative correction turns to zero 
at a certain value of r as the result of the cancellation between the quark condensate 
contribution (~ const) and that of the mixed condensate (~ r). This cancellation of 
course cannot be taken seriously and only indicates that the OPE breaks down since the 
hierarchy of contributions is lost. We conclude that the classical QCD sum rule is only 
valid up to light-cone distances of order 1 — 3 GeV~^, depending on the value of the Borel 
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Figure 6: Perturbative (solid curves) and nonperturbative (long dashes) contributions to the 
B-meson distribution amplitude ^+{t, = 1 GeV) calculated from the sum rule to the NLO 
accuracy. In addition, the nonlocal condensate models pSj) . H36|) of resummed nonperturbative 
contributions to the sum rule, cf. (|37|) . are shown by short dashes. The continuum threshold 
is chosen to be wo=l GeV and two values of the Borel parameter are used: M = 0.3 GeV (left 
panel) and M = 0.6 GeV (right panel). The value of the decay constant -F(^) appearing on the 
l.h.s. of is substituted by the corresponding sum rule 
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Figure 7: Same as in Fig. 6, with a different value of the continuum threshold = 0.8 GeV. 



parameter. A rough estimate for the nonperturbative contribution to in the strict 
OPE-based approach is, therefore, given by the integral over the region of small r where 
the correction is still positive, that is up to the crossing point with the zero axis. In 
order to get an estimate of a possible nonperturbative contribution from large distances 
we use the concept of a nonlocal quark condensate introduced in jSI] and later used rather 
extensively in QCD sum rule calculations of the distribution amplitudes of light mesons 
by the Dubna group [221 EBj- The same approach was taken up in 14j. 

The nonlocal quark condensate presents a model for a partial resummation of the OPE 
to all orders in terms of the vacuum expectation value of the single nonlocal operator 

/■oo 

(0|g(x)[x,0]g(0)|0) = (gg) / du e^^Vl^) • (33) 

Jo 

The first two moments of f{v) are fixed by the OPE: 

dvf{v) = l, duufiu) = -ml, (34) 
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and in addition one requires that the correlation function ()33p decreases exponentially at 
large Euclidian separations —>■ — oo. The two simplest choices are 

Model I : f{u) = 6{iy - ml/A) ^ (35) 

corresponding to the Gaussian large-distance behavior ~ exp[— |x^|mQ/16] and 

Model II: f(u) = ——-y-Pe''/'', P = 3 + — |2H1 (36) 



corresponding to (0|g(a;)g(0)|0) ~ exp[— Ay^x^]. Here A is a parameter with physical 
meaning of the vacuum quark correlation length. In this work we take A = 400 MeV as 
a representative number, cf [321 El] • We will see that sensitivity of the sum rules to the 
shape of f{v) is in fact small; the major shortcoming of this approach is rather that other 
condensates (e.g. the nonlocal quark-antiquark-gluon condensate) are not included and 
there exists no parameter that would justify such a truncation. 

Nonlocality of the quark condensate is easy to implement within our sum rules and it 
amounts to a simple substitution in Eq. (|3^ (cf. i l4j ) 

--{qq) S{t, M, /.)|i + 0{as) - -j^ (1 + 2rM) 

----{qq)S{r,M,^,) dz. /(z.) e^^^+^^^V^^^^) . (37) 

Note that the mixed condensate contribution is now included as a part of the nonlocal 
condensate and we neglect the ttg-correction to the (local) quark condensate (but retain 
the Sudakov exponent). The results are shown in Fig. 6 and Fig. 7 by short dashes; the 
lower and the higher of the curves correspond to the choice in ()35|) and ()36|). respectively. 
The corresponding results for A^^ are, for = 1 GeV: 



UJQ 



1 GeV , M = 0.6 GeV : A;b^ = 1.23 + { 0.60 = 1.95 ± 0.23 GeV"^ , 



LUo = 1 GeV , M = 0.3 GeV : Ao^ = 1.32 + I 0.54 = 2.03 ± 0.29 GeV"^ , 



cuo = 0.8 GeV , M = 0.6 GeV : Ao^ = 1.36 + { 0.84 = 2.36 ± 0.33 GeV"^ , 



ujo = 0.8 GeV , M = 0.3 GeV : A;b' = 1.39 + < 0.64 = 2.24 ± 0.35 GeV"\ (38) 

[ 1.05 

where the first number gives the perturbative contribution (the difference to (j22j) is due 
to the different value used for -F(/i)) and the three numbers under the brace correspond 
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to three different estimates for tlie nonperturbative contribution: 1) quark and mixed 
condensate contribution to ()32|) restricted to tlie positivity region, 2) nonlocal conden- 
sate model I (j35|) and 3) nonlocal condensate model II (j36|l . The first (upper) number 
should be considered as an estimate of the nonperturbative correction from below while 
the difference between the two lower ones characterizes the uncertainty in the choice of the 
parametrization of the nonlocal condensate. We take the average between the two models 
as our central value, and one half of the difference between this central value and the first 
(upper) number, coming from local OPE, as an estimate of the overall uncertainty of the 
result. In other words, we ascribe 50% uncertainty to the extrapolation of the nonper- 
turbative contribution to large distances as suggested by the nonlocal condensate model, 
which is rather conservative, cf. [22] • From the numbers in Eq. (jHHj) we obtain the final 
result 

A;b^(/x= 1 GeV) = 2.15 ±0.5 GeV"^ (39) 

or 

Xsifi = 1 GeV) = 460 ± 110 MeV . (40) 

Our value of is somewhat larger than the number accepted in jSHl E] = 0.35 ± 
0.15 GeV, although consistent with it within errors. It is also consistent with the rough 
estimate A^ ~ 0.6 GeV derived in [15J. 

As follows from (jSj) and ©, the scale dependence of Xb also involves the first logarith- 
mic moment of the distribution amplitude fHl 1^ 



asCp /i 

1 H m — 

2tt fiQ 



A.^(/.o)-^ln^ rf lnf^0.(fc,/.o), (41) 
71 fJ^o Jo k k 



where (ttsCir/vr) ln(yu//io) < 1. We define 

dk u 

— In ^Mk,t^) = Mt^) J ^ dr ln(r/ieT^)(^+(r,/i) (42) 

and calculate o"b(1 GeV) from the QCD sum rule (j32|) repeating the same procedure as 
explained above for A^. Without going into details we simply quote the final result 

(TB(/i = 1 GeV) = 1.4 ±0.4. (43) 

Note that crsifJ') defines the average value of In/i/fc in the integral for the first inverse 
moment A^^, so that the number in ()43|) implies that main contribution to A^^ comes 
from momenta ~ 250 MeV. With this value for as, the two contributions 0{as) in (j4ip 
tend to cancel each other to a large extent, so that the remaining scale dependence of A^^ 
is weak. 

A simple model of the distribution amplitude (f)+{k, fi) with given values of the param- 
eters Xb and as and correct asymptotic behavior can be chosen as 

(/.+ (A;,^ = 1 GeV) - ^ 



71 k^ + 1 



k^ + 1 7r2 



(44) 



{k in units of GeV). Using the values of Ag and ctb in (jHU)) . one obtains the distribution 
shown in Fig. 8 by the solid curve. For comparison, we also show in this plot a typical 
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12 3 4 

k, GeV 

Figure 8: A QCD model for the B-meson distribution amplitude (|44() (solid curve) compared 
with the perturbative sum rule prediction with M=0.45 GeV, ojq = \ GeV (dashed curve). 

distribution amplitude obtained from QCD sum rules in perturbation theory (fT^ in Sec. 3, 
cf. Fig. 3. The effect of nonperturbative corrections is to shift the distribution towards 
softer momenta, which is natural. One minor drawback of such a parametrization is that 
the set of parameters and cr^ is not closed under renormalization. In view of a very 
limited range of scales that are interesting for B-decay phenomenology this seems to be 
not a problem, however. 

To summarize, in this paper we have derived QCD sum rules for the B-meson distribu- 
tion amplitude (Q) and, in particular, obtained an estimate of its first inverse moment A^^ 
and the parameter as that characterizes the shape of the distribution, see (jlBI). 
A simple model is suggested (j44j) that incorporates all existing constraints. We believe 
that our estimates are interesting for the studies of the heavy quark limit in exclusive 
B-decays and can be used in a broad context. Concrete applications go beyond the task 
of this work. 
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